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O ' Abstract 
O 

o. 

An astrophysically realistic model of wave dynamics in black-hole spacetimes 

g. 

must involve a non-spherical background geometry with angular momentum. 



(N 



We consider the evolution of gravitational (and electromagnetic) perturbations 



black hole. Moreover, our results turn over the traditional belief (which has 



. in rotating Kerr spacetimes. We show that a rotating Kerr black hole becomes 

ON . "bald" slower than the corresponding spherically-symmetric Schwarzschild 

t> 
o 

o\ . 

Q\ . been widely accepted during the last three decades) that the late-time tail of 

*0 ■ 

Q^. gravitational collapse is universal. In particular, we show that different fields 

bJT). have different decaying rates. Our results are also of importance both to the 



study of the no-hair conjecture and the mass-inflation scenario (stability of 
Cauchy horizons). 

The radiative tail of gravitational collapse decays with time leaving behind a Kerr- 
Newman black hole characterized solely by the black-hole mass, charge, and angular- 
momentum. This is the essence of the no-hair conjecture, introduced by Wheeler in the 
early 1970s 0. 

Price was the first to analyze the mechanism by which the spacetime outside a 
(nearly spherical) star divests itself of all radiative multipole moments, and leaves behind a 
Schwarzschild black hole; it was demonstrated that all radiative perturbations decay asymp- 
totically as an inverse power of time, the power indices equal 2/ + 3 (in absolute value), where 
I is the multipole order of the perturbation. This late-time decay of radiative fields is often 
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referred to as their "power-law tail" . Physically, these inverse power-law tails are associated 
with the backscattering of waves off the effective curvature potential at asymptotically far 
regions 

The analysis of Price has been extended by many authors. We shall not attempt to review 
the numerous works which address the problem of the late-time evolution of gravitational 
collapse. For a partial list of references, see e.g., Q. These earlier analyses were restricted, 
however, to spherically symmetric backgrounds. It is well-known that realistic stellar objects 
generally rotate about their axis, and are therefore not spherical. Thus, the nature of the 
physical process of stellar core collapse to form a black hole is essentially non-spheric. An 
astrophysically realistic model must therefore take into account the angular momentum of 
the background geometry. 

The corresponding problem of wave dynamics in realistic rotating Kerr spacetimes is 
much more complicated due to the lack of spherical symmetry. A first progress has been 
achieved only recently [§]-§]. Detailed analyses for the simplified toy-model of a test scalar 
field in the Kerr background have been given recently in PJIU|. 

Obviously, the most interesting situation from a physical point of view is the dynamics 
of gravitational perturbations in rotating Kerr spacetimes. This is the subject of this Letter, 
in which we present our main results for this fascinating problem. Full details of the analysis 
are given elsewhere H|. 

The dynamics of massless perturbations outside a realistic rotating Kerr black hole is 
governed by Teukolsky's master equation 



(r 2 + a 2 ) 2 



A 

d 



a 2 sin 2 9 



d 2 ib AMar d 2 ib 
+ I n J + 



dt 2 



A dtdcp 



-A~ s — A s — 
dr \ dr 

M(r 2 - a 2 ) 



1 d ( . df 

sine 09 [ s de 



-2s 



a i 
A ~ sin 2 9 
air — m 



d 2 i) 

dip 2 
i cos 6 

+ 



sin 2 6 



dip 



r — ia cos 9 



^ + (s 2 cot 2 9-s)?P = 
at 



(1) 



where M and a are the mass and angular-momentum per unit-mass of the black hole, and 
A = r 2 — 2Mr + a 2 . (We use gravitational units in which G — c — 1). The parameter s is 



called the spin-weight of the field. For gravitational perturbations s = ±2 (for electromag- 
netic perturbations s = ±1). The field quantities ip which satisfy Teukolsky's equation are 
given in [|12[ . 

oo 

Resolving the field in the form ip = A~ s ' 2 (r 2 + a 2 ) -1 ' 2 J2 ^> rn e vmip (where m is the 

m=— oo 

azimuthal number), one obtains a wave-equation for each value of m 
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where the tortoise radial coordinate y is defined by dy = [(r 2 + a 2 )/A]dr. (We suppress the 
index m). The coefficients Bi = Bi(r, 9) are given by 
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[The explicit expression of B 3 (r,9) is not important for the analysis]. 
The time-evolution of a wave-field described by Eq. (|2|) is given by 
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for t > 0, where z stands for (y, 9). The (retarded) Green's function G(z, z'\ t) is defined by 
DG(z, z'\ t) = 5{t)5(y - y')5{9 - 0')/27rsin0 with G = for t < 0. We express the Green's 
function in terms of the the Fourier transform Gi(y, y'; w) 
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where c is some positive constant and l = max(\m\, \s\). The functions s Sl n (9,aw) are the 
spin-weighted spheroidal harmonics which are solutions to the angular equation Jl2| 
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The Fourier transform is analytic in the upper half w-plane and it satisfies the equation 
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wherein = (r 2 +a 2 )w—am, X = A+a 2 w 2 —2amw, and H = s(r—M)/(r 2 +a 2 )+rA/(r 2 +a 2 ) 2 . 

Define two auxiliary functions and as solutions to the homogeneous equation 
D(w)i&i=D(w)i&2 — 0, with the physical boundary conditions of purely ingoing waves cross- 
ing the event horizon, and purely outgoing waves at spatial infinity, respectively. In terms 
of \&i and \^2, and henceforth assuming y' < y, Gi(y,y';w) = —^fi(y',w)^f2(y,w)/W(w), 
where we have used the Wronskian relation W(w) = W^i, ^2) = ^i^2, y — ^2^i,y 

It is well known that the late-time behaviour of massless perturbations fields is deter- 
mined by the backscattering from asymptotically far regions Thus, the late-time 
behaviour is dominated by the /ou^frequencies contribution to the Green's function, for only 
low frequencies will be backscattered by the small effective curvature potential (at r ^> M) . 
Therefore, a small-w approximation [or equivalently, a large-r approximation of Eq. (^|)] is 
sufficient in order to study the asymptotic late-time behaviour of the fields |0 . With this ap- 
proximation, the two basic solutions required in order to build the Fourier transform are ^1 = 
r l+1 e iwT M{l+s + l-2iwM, 21+2, -2iwr) and f 2 = r l+1 e iwr U(l+s+l-2iwM, 21+2, -2iwr), 
where M(a, b, z) and U(a, b, z) are the two standard solutions to the confluent hypergeomet- 
ric equation 0. Then W(*i,* 2 ) = i{-l) l+1 {2l + \)\{2w)-^ 2l+ ^ / {I + s)\. 



In order to calculate G(z,z';i) using Eq. (|6]), one may close the contour of integration 
into the lower half of the complex frequency plane. Then, one identifies three distinct contri- 
butions to G{z, z'\ t) |L5| : Prompt contribution, quasinormal modes, and tail contribution. 
The late-time tail is associated with the existence of a branch cut (in ^2) in the complex fre- 
quency plane |HJ (usually placed along the negative imaginary w— axis). A little arithmetic 
leads to ITBl 
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(2/+1)! W(w) 
Taking cognizance of Eq. (|6]) we obtain 

c - l M(-iy2^(l + s)\(l-s)\ 

/ ^(y, w)y 1 (y',w) s S l (e,aw) s S l (6',aw)w 2l+2 e- lwt dw . (10) 
Jo 

The angular equation (0) is amenable to a perturbation treatment for small aw; we write 
it in the form (L° + L 1 ) s S] n = —gA^sS™, where ^(6, aw) = (aw) 2 cos 2 6 — 2awscos6 [and 
L°(6 ) ) is the w-independent part of Eq. (0)], and we use the spin-weighted spherical functions 
gY™ as a representation. They satisfy L° s Yi = — s Ai S Yi with S A| ' ) = (i — s) (I + s + 1). For 
small aw a standard perturbation theory (see, for example, |T7[]) yields [|J 



a Si{6 } aw) = £ ^(a^i'^ir,^) , (11) 

fc=Z 

where, to leading order in aw, the coefficients Cik(aw) are w— independent ||. Equation 
(|TT|) implies that the black-hole rotation mixes different spin-weighted spherical harmonics. 

The time-evolution of the fields is given by Eq. (^). Therefore, in order to eluci- 
date the coupling between different modes we should evaluate the integrals {slm\skm), 
(slm\ sin 2 8\skm), and (slm\ cos 8\skm), where (slm\F(6)\skm) = j s Y[ Tn *F(6) s Y™dQ [see 
Eqs. (^|) and (|j) for the definition of the Bi(r,8) coefficients]. The spin- weighted spher- 
ical harmonics are related to the rotation matrix elements of quantum mechanics [JTSj . 



Hence, standard formulae are available for integrating the product of three such functions 
(these are given in terms of the Clebsch-Gordan coefficients H). In particular, the inte- 
gral (s/0| sin 2 9\sk0) vanishes unless I = k, k ± 2, while the integral (sl0\ cos#|s£;0) vanishes 
unless I — k ± 1. For non-axially symmetric [m ^ 0) modes, (slm\ sin 2 6\skm) ^ for 
I = k, k ± 1, k ± 2, and (slm\ cos9\skm) ^ for I = k, k ± 1 (all other matrix elements 
vanish) . 



Asymptotic behaviour at timelike infinity. As already explained, the late-time behaviour 
of the fields should follow from the Zoiy-frequency contribution to the Green's function. 
Actually, it is easy to verify that the effective contribution to the integral in Eq. fllQ|) should 
come from \w\=0(l/t). Thus, in order to obtain the asymptotic behaviour of the fields at 
timelike infinity (where y,y' <C t) we may use the \w\r <C 1 asymptotic limit of ^i(r,w), 
which is given by \E'i(r, w) ~ r l+1 [Q. 

Substituting this in Eq. (|T0|) , and using the representation Eq. ([TT]) for the spin- weighted 
spheroidal wave functions s Si [togather with the above cited properties of the integrals 
(slm\ sin 2 6\skm) and {slm\ cos9\skm)}, we find that the asymptotic late-time behaviour of 
the I mode (where / > Iq) is dominated by the following effective Green's function: 

Gf(z, z'; t) = MF^yyy^KY^^ey^H-^ 1 ^ , (12) 

where F\ = Fi(l, l , m, s) = (-l)( w »+ 2s + 2 )/ 2 2 2 ^Z + Z + 2^^ 

We emphasize that the power indices I + Iq + 3 found here for rotating Kerr spacetimes are 
smaller than the corresponding power indices (the well known 21+3) in spherically symmetric 
Schwarzschild spacetimes. (There is an equality only for the I = Iq mode). This implies a 
slower decay of perturbations in rotating Kerr spacetimes. 

Asymptotic behaviour at future null infinity. It is easy to verify that for this case the 
effective frequencies contributing to the integral in Eq. ([K]) are of order 0(l/u). Thus, 
for y — y' <C t <C 2y — y' one may use the \w\y' <C 1 limit for ^i(y',w) and the M <^ 
|u>| _1 y (Imw < 0) asymptotic limit of ^i(y,w), which is given by ^i(y,w) ~ e twy (2l + 

1 jj e -i ff (i+H-l)/2( 2t u)-(J+»+l)j,--/(/ - S )\ 0. 

Substituting this in Eq. (|T0|), and using the representation Eq. (|TTD for the spin- weighted 
spheroidal wave functions, we find that the behaviour of the / mode (where / > Iq) at 
the asymptotic region of null infinity scri + is dominated by the following effective Green's 
function: 

Gf(z, z'; t) = J2 MF 2 y' k+1 v- s s Y l (9) s Y k *(9')a l - k u-^- s+ ^ , (13) 

k=lo 
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where F 2 = F 2 (l, k, m, s) = (-l)(«+M-2.+2)/2 2 k(/. + s yy _ s + i)!C , fc ,/ 7r (2A; + 1)!. 

Asymptotic behaviour at the black-hole outer horizon. The asymptotic solution to the 
homogeneous equation D(w)^i(y,w) = at the black-hole outer horizon H + (y — ► — oo) 
is = C( w )A- s / 2 e- i ^- m,u +^ Q, where w+ = a/2Mr + [r+ = M + (M 2 - a 2 ) 1 / 2 

is the location of the black-hole outer horizon]. In addition, we use ^\{y', w) ~ y +1 . 
Regularity of the solution requires C to be an analytic function of w. We thus expand 
C(w) = C + CiW + ■ • • for small w (as already explained, the late-time behaviour of the 
field is dominated by the /ou>-frequency contribution to the Green's function). Substituting 
this in Eq. fllUP , and using the representation Eq. (|TT| ) for the spin- weighted spheroidal 



wave functions, we find that the asymptotic behaviour of the I mode (where I > Iq) at the 
black-hole outer horizon H + is dominated by the following effective Green's function: 

Gf(z, z'; t) = s r|Mf 1 A- s / 2 2/' i ° +1 ^(^^(^> | -'°e ,ra '"+ ! ');- (l+io+3+6) , (14) 

where S T[ are constants, and 6 = generically, except for the unique case m = with s > 0, 
in which b = 1 @. 

Pure initial pulse. So far we have assumed that the initial pulse consists of all the allowed 
(I > Iq) modes. If, on the other hand, the initial angular distribution is characterized by 
a pure spin- weighted spherical harmonic function sY™, then the asymptotic late-time tails 
are dominated by modes which, in general, have an angular distribution different from the 
original one (a full analysis of this case is given in We find that the field's behaviour 
at the asymptotic regions of timelike infinity i + and at the black-hole outer horizon H + is 
dominated again by the lowest allowed mode (i.e., I = Iq). The damping exponents are (in 
absolute value) I* + Iq + 3 — q and I* + Iq + 3 — q + b, respectively, where q = min(l* — Iq, 2). 

On the other hand, the behaviour of gravitational (and electromagnetic) perturbations at 
the asymptotic region of null infinity scri + is dominated by the I = Iq mode if Iq < I* < Iq + 2 
and by the Iq < I < I* — 2 modes otherwise The corresponding damping exponents are 
Iq — s + 2 and I* — s, respectively. 

Summary and physical implications. We have analyzed the dynamics of gravitational 



(physically, the most interesting case) and electromagnetic fields in realistic rotating black- 
hole spacetimes. The main results and their physical implications are as follows: 

(1) We have shown that the late-time evolution of realistic rotating gravitational collapse 
is characterized by inverse power-law decaying tails at the three asymptotic regions: timelike 
infinity i + , future null infinity scri + , and the black-hole outer-horizon H + (where the power- 
law behaviour is multiplied by an oscillatory term, caused by the dragging of reference 
frames at the event horizon). The relaxation of the fields is in accordance with the no-hair 
conjecture 0. This Letter reveals the dynamical physical mechanism behind this conjecture 
in the context of rotating gravitational collapse. 

(2) The unique and important feature of rotating gravitational collapse is the active 
coupling between modes of different I (but the same m). Physically, this phenomena is 
caused by the dragging of reference frames, due to the black-hole (or star's) rotation (this 
phenomena is absent in the non-rotating a = case). As a consequence, the late-time 
evolution of realistic rotating gravitational collapse has an angular distribution which is 
generically different from the original angular distribution (in the initial pulse). 

(3) The power indices at a fixed radius are found to be / + l + 3. These damping expo- 
nents are generically smaller than the corresponding power indices in spherically symmetric 
spacetimes. This implies a slower decay of perturbations in rotating Kerr spacetimes. Stated 
in a more pictorial way, a rotating Kerr black hole generically becomes "bald" slower than 
a spherically-symmetric Schwarzschild black hole. 

(4) It has been widely accepted that the late-time tail of gravitational collapse is uni- 
versal in the sense that it is independent of the type of the massless field considered (e.g., 
scalar, neutrino, electromagnetic, and gravitational). This belief was based on spherically 
symmetric analyses. Our analysis, however, turnover this point of view. In particular, the 
power indices / + / + 3 at a fixed radius which are found in this Letter are generically different 
from those obtained in the scalar field toy-model ||[10| / + \m\ + p + 3 (where p — if I — \m\ 
is even, and p = 1 otherwise). 

We have shown that different types of fields have different decaying-rates. This is a 
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rather surprising conclusion, which has been overlooked in the last three decades ! It should 
be stressed, therefore, that the results obtained from the scalar field toy-model [0,|1OJ are 
actually not applicable for the physically interesting case of higher-spin perturbations (i.e., 
gravitational and electromagnetic fields). 

(5) Our results should have important implications for the mass- inflation scenario and the 
stability of Cauchy horizons (see e.g., |20|j2l[1 and references therein). In particular, the late- 



time tails found in this Letter should be used as initial data for perturbations propagating 
inside the (rotating) black hole. 
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